Band structure of Charge Ordered Doped Antiferromagnets. 
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We study the distribution of electronic spectral weight in a doped antiferromagnet with various 
types of charge order and compare to angle resolved photoemission experiments on lightly doped 
La2-iSr2,Cu04 (LSCO) and electron doped Nd2-a;Cea;Cu04±i. Calculations on in-phase stripe and 
bubble phases for the electron doped system are both in good agreement with experiment including 
in particular the existence of in-gap spectral weight. In addition we find that for in-phase stripes, 
in contrast to anti-phase stripes, the chemical potential is likely to move with doping. For the 
hole doped system we find that "staircase" stripes which are globally diagonal but locally vertical 
or horizontal can reproduce the photoemission data whereas pure diagonal stripes cannot. We 
also calculate the magnetic structure factors of such staircase stripes and find that as the stripe 
separation is decreased with increased doping these evolve from diagonal to vertical separated by a 
coexistence region. The results suggest that the transition from horizontal to diagonal stripes seen 
in neutron scattering on underdoped LSCO may be a crossover between a regime where the typical 
length of straight stripe segments is longer than the inter-stripe spacing to one where it is shorter 
and that locally the stripes are always aligned with the Cu-0 bonds. 



I. INTRODUCTION 

In several families of superconducting cuprates there 
is evidence for stripes which are regularly spaced quasi- 
one dimensional structures where the doped charge 
assembles.^ At the same time there is a well developed 
theory of high-temperature superconductivity in a sys- 
tem of weakly coupled Hubbard ladders uSi^ It is thus 
quite natural to speculate that such a theory is in fact 
realized in the cuprates. One obvious objection to the 
stripe scenario of superconductivity is the lack of con- 
vincing evidence for the existence of stripes in several 
materials given that the strongest evidence is found in the 
relatively low-Tc La2_2:Sra;Cu04 (LSCO) family.^ From 
a theory point of view the lack of direct evidence is not 
immediately discouraging because stripe order is bad for 
superconductivity whereas more elusive dynamic charge 
stripe correlations are good.--^ However, even in a sys- 
tem such as LSCO where stripes are well established 
there are still issues about their fundamental implica- 
tions on the electron dynamics. In particular, it seems 
the distribution of low-energy spectral weight in k-space 
as probed by angle-resolved photoemission spectroscopy 
(ARPES) does not show any clear evidence of the sug- 
gested quasi-one-dimensional nature of the electronic 
states, looking instead like the Fermi surface of a fully 
two-dimensional systemn^ On the other hand ARPES also 
provides clear evidence of exotic physics, perhaps elec- 
tron fractionalization^, with spectral functions that are 
typically very broad in energy and not consistent with a 
Fermi liquid based quasi-particle descriptioniiS 

A simple model for studying the distribution of sin- 
gle electron spectral weight in a charge ordered doped 
antiferromagnet was introduced by Salkola et alJ^^ in 
which all the complicated physics which is responsible for 
the stripe formation and integrity is replaced by a hand 



picked potential which emulates the local environment of 
electrons in a stripe ordered antiferromagnet. The poten- 
tial is simply a staggered field representing the local anti- 
ferromagnetism together with anti-phase domain walls of 
suppressed field strength representing stripes. The spac- 
ings of the stripes are chosen in such a way that the 
model by construction will reproduce the diffraction re- 
sponse of a stripe ordered system. The stripe placements 
are static but may be chosen in an irregular fashion sim- 
ulating quenched disorder or dynamic stripes which are 
fluctuating slowly compared to time scales of the local 
electron dynamics. Using this model for a disordered 
array of quarter filled "bond-aligned" ■'^^ stripes it was 
found that the low energy spectral weight forms a two- 
dimensional Fermi surface with, in particular, spectral 
weight in the "nodal regions", near (7r/2,7r/2), which 
naively corresponds to propagation diagonally with re- 
spect to the stripe direction. A more detailed study of 
this model was carried out in Ref. where it was real- 
ized that even a single stripe with one-dimensional states 
localized transverse to the stripe captures the qualita- 
tive features of the low-energy spectral weight as seen in 
ARPES in the under and optimally doped LSCO with 
the full spectral weight of the stripe lying roughly on 
the diamond Fermi surface of a half-filled nearest neigh- 
bor tight-binding model. Here also issues of interactions 
along a stripe was addressed and it was found that the 
distinct non- Fermi liquid properties of an interacting one- 
dimensional electron gas may be displayed also in direc- 
tions not aligned with the stripes. 

In the present paper we extend the work of Ref. ^ and 
IT^ to look at a broader range of ordered structures. This 
is motivated by a search for a more comprehensive under- 
standing of charge order in a doped antiferromagnet but 
also by direct or indirect observations of structures which 
are not consistent with bond-aligned anti-phase stripes. 
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In the very lightly doped, non-superconducting, phase 
of LSCO there are diffraction patterns consistent with 
diagonal stripesi^^ instead of the bond-aligned stripes 
seen at higher doping. The transition between bond- 
aligned and diagonal stripes coincides with the supercon- 
ductor to insulator transition at doping x w 0.055, sug- 
gesting a close link between the two properties. Within 
the non-interacting model considered here we cannot di- 
rectly address the connection between superconductivity 
and stripe orientation but by comparing the predicted 
spectral weight of diagonal stripes with that observed 
in ARPE^S, we gain some insight into the nature of 
such diagonal stripes. Recently, there has also been indi- 
rect evidence both from NMR^^ and thermal conductiv- 
ity measurements^^ of an inhomogeneous charge struc- 
ture in the electron-doped cuprates (Pr,La)2-3;Cea;Cu04 
(PLCCO). However, in electron doped materials there 
has been no evidence for incommensurate magnetismiS 
effectively ruling out the possibility of anti-phase stripes. 
For this reason we have studied in-phase structures and 
compared to ARPES results on the Nd2-2;Cej;Cu04±5 
(NCCO) family^ 

Our results can be summarized as follows. As exempli- 
fied in Fig. [21we find that the spectral weight of localized 
states is centered on the Brillouin zone (BZ) diagonals, 
given by cos(/i:a;)-|-cos(fcj^) — 0, independently of the mag- 
nitude of second and third nearest neighbor hopping t' 
and t" as well as the shape and form of the charge "im- 
purity" . This shows why in general one may not expect 
any dramatic signature of stripes or other inhomogeneous 
electronic structure on the k-space distribution of low en- 
ergy spectral weight, as this is where the Fermi surface 
is located within a nearly half-filled tight-binding model 
which is dominated by nearest neighbor hopping. How- 
ever, as shown in Fig. [7| El and El a defining feature 
of the stripe states is that they have an energy which 
is within the Mott gap of the undoped system. Such 
in-gap states appear to be a common feature of the evo- 
lution of the band structures which doping as measured 
in ARPES and which we believe is a strong indication of 
an inhomogeneous charge distributionn^i^ 

We compare the dispersion and distribution of spectral 
weight of diagonal and bond-aligned stripes and find that 
pure diagonal stripes cannot reproduce the characteris- 
tic "Fermi arc" centered around the nodal region which is 
seen in ARPES on lightly doped LSCO. Instead we find 
that the spectral weight of a hole doped diagonal stripe 
is concentrated to the "anti-nodal" BZ regions around 
(tt, 0) with very little weight in the nodal region. In ad- 
dition, the band width of states on a diagonal stripe is 
expected to be roughly 2\t'\ < 0.2eV which is inconsistent 
with the ARPES data where the band width of the in- 
gap states is of the order of leV (assuming half is seen). 
However, it turns out that the spectral distribution and 
band width of bond-aligned stripes is qualitatively con- 
sistent with the ARPES data. For this reason we suggest 
that the diagonal stripe phase consists of stripes which 
are globally diagonal but locally bond-aligned, a carica- 



ture of which are the "staircase stripes" shown in Fig. 01 
In Fig. Ois shown the low energy spectral weight which 
is concentrated around the nodal region in a hole doped 
array of staircase stripes. 

An interesting aspect of these staircase stripes is the 
magnetic structure factors which depend on the ratio 
between the length of the bond-aligned segments, the 
"step" length I, and the distance between neighboring 
stripes d. As shown in Fig. Owe find that the correspond- 
ing Bragg peaks can be classified in three main regimes. 
For / « d there are two peaks corresponding to diagonal 
stripes along the x + y direction at (vr ± ^diag, tt T ^dmg), 
for I « 2d the diagonal peaks coexist with four peaks at 
(tt, 7r±(5coi) and {nzLScoi, tt) corresponding to the response 
expected from bond-aligned stripes along the x and y di- 
rections, while for I ^ d there are only the four bond- 
aligned peaks but they are now shifted away slightly 
from the square lattice axes. The qualitative features 
are very similar to what is found from neutron scattering 
in LSCO as a function of dopingi^, although the rela- 
tive incommensurability (5diag/<^coi is not quite accurately 
reproduced in the coexistence regime. Nevertheless, this 
suggests a scenario in which the stripes in the orthorhom- 
bic phase of LSCO are always locally bond-aligned but 
in some sense globally diagonal with a crossover as the 
stripe spacing is decreased with increased doping and 
not a first order transition as suggested by the neutron 
scattering data. Corroborating such a crossover scenario 
is the fact that the ARPES spectra evolve smoothly 
through the diagonal to bond-aligned stripe transition 
and that the Fermi velocity in the nodal direction is 
roughly independent of doping, indicating that if the low- 
energy spectral weight is stripe related the local char- 
acter of the stripes does not change dramatically with 
doping.^*' In addition, the hole mobility at moderate tem- 
peratures changes by only a factor of three from very light 
{x = 0.01) to optimal {x = 0.17) doping which is very 
naturally understood within a stripe model in which the 
local character of the stripes is roughly independent of 
doping»2^ 

In Fig. El we show the the band structure of a sys- 
tem with disordered in-phase stripes and in Fig. Illl the 
corresponding electron doped low energy spectral weight 
which is in qualitative agreement with the ARPES data. 
The most interesting part of these results is the evolution 
of spectral weight in the nodal region, where at light dop- 
ing (4%) there is in-gap spectral weight which at higher 
doping (10%) broadens as a consequence of the shorter 
inter-stripe distance and reaches the Fermi surface. The 
evolution of the Fermi surface with doping can be repro- 
duced in mean-field theory of Hubbard model with longer 
range hopping by allowing for a doping dependent inter- 
action U.'^'^ The difference within a stripe model is that 
the low-energy states are dynamically one-dimensional, 
being localized transverse to the stripes (Fig. I13() and 
the existence of in-gap spectral weight. In addition we 
find that for in-phase stripes, in contrast to anti-phase 
stripes, the chemical potential is likely to move into the 
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upper Hubbard band with electron doping because the 
in-phase stripe states he close to the upper and lower 
Hubbard bands (Fig. 

Finally, we present similar results from a calculation 
on a "bubble" phase where the doped charge is con- 
fined to small "O-dimensional" droplets instead of the 
1-dimensional stripes. Bubbles would arise naturally in- 
stead of stripes in a i — J model with long-range Coulomb 
repulsion in the limit t <C J because of the lower mag- 
netic energy. Thus in the electron doped materials, which 
appear to have "stronger" antiferromagnetism than the 
hole doped materials, one may speculate that bubbles are 
favored over stripes. As far as the distribution of spectral 
weight is concerned (Fig. [T5land [THjl there is little qual- 
itative difference between bubbles and stripes. For the 
bubbles we find that the nodal spectral weight broadens 
as the size of the bubbles increase with doping in analogy 
with the increasing density of stripes. 



II. THE MODEL 

We will consider a tight-binding model on a square lat- 
tice with first, second and third nearest neighbor hopping 
together with a static potential which represents stripes 
or other charge structures. The Hamiltonian reads 

H = - tJ2 (4,.CrV + H.C.) 

{rr'){T 

- t' {cl.aCr'a+H-G.) 

{rr')' (7 

- t" (4,c,v+H.C.) 

{rr')" a 

+ a(-ir+^F(a:,y)4^,,c,,,.., (1) 

where Cr^a is the electron destruction operator at site 
r = {x,y) and with spin cr = ±. The hopping is given in 
a standard fashion where {rr') indicates nearest neigh- 
bors, {rr'Y next-nearest neighbors, and {rr')" next-next 
nearest neighbors. In what follows we will use energy 
units such that i = 1 and we take t' jt < and t" jt > 0. 
Unless stated otherwise we will be using m = t = 1 . The 
physical intuition for the field m is that it is the energy 
cost of moving a hole from a stripe into the antiferro- 
magnetic background and it is thus expected to be of the 
order of the spin exchange J. We also use units such that 
the lattice constant is equal to one and h = 1. 

The potential m{—iy^~^yV{x,y) describes the collec- 
tive field which defines the stripe order. We will use the 
simplest possible form and take V{x, y) =1,0 or -1, where 
V{x,y) —1 or -1 represent antiferromagnetic regions re- 
lated by a TT phase shift and V{x^ y) = are the locations 
of the stripes. The case V{x, y) = 1 for all x and y corre- 
sponds to the standard mean-field result of the Hubbard 
model at half-filling^'* giving an upper and lower Hubbard 
band separated by the Mott-Hubbard gap. Introducing 



regions where V{x, y) = will in general give rise to 
localized "impurity" states within the gap. 

The object which was studied in detail in Ref 0| is the 
site-centered anti-phase stripe given by 



stripe 



{x,y) 




(2) 



and displayed graphically in Fig. Here we also show 
two other charge structures, diagonal stripes and bub- 
bles, which we will consider in more detail subsequently. 
In a real system we want to consider ordered or disor- 
dered arrays structures such that at some finite doping 
the corresponding states are partly occupied. 



(a) 



FIG. 1: Graphic representation of the potential V{x,y) for 
various charge structures. Here black corresponds to = 0, 
gray to V — 1, white to V — —1 and each square represents 
a site {x,y) on the lattice, (a) is a bond-aligned anti-phase 
stripe, (b) is an anti-phase diagonal stripe and (c) a bubble. 



A. Spectral weight of localized states 

Here we will show that a state localized on an impurity 
where the staggered field is zero has its spectral weight 
centered within a region m of the cos(fc^) + cos(fcj,) = 
diamond, regardless of the values of the longer range hop- 
ping t' and t" and the geometry of the impurity. This 
result is really a trivial consequence of the fact that the 
nearest neighbor hopping t connects the two sublattices 
of the staggered field while t' and t" does not, or equiva- 
lently the symmetry or not of the dispersion with respect 
to a shift of the momenta by scattering vector (tt, tt) of 
the staggered field. 

We define a potential for an arbitrary impurity 

V{x, y) — m{~\)^^y , {x, y} ^ impurity 
V{x,y)^ 0, e impurity, (3) 

and write the tight-binding dispersions 



= — 2^(cos(fca;) + cos{ky)) 



(4) 



-At' cos{kx) cos{ky) — 2t" {cos{2kx) + cos{2ky)) , 



With this we solve for the eigenstates in the bulk 



(5) 
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with energy 

E-^ = e\±^(e^^^ + m\ (6) 

and ratio of coefficients 

= {-el ± ^(£0)2 + m2)/m . (7) 

The locahzed state can be expanded in terms of the 
complete set of states ijjk which are the solutions in the 
bulk. Clearly any state which is localized at the impurity 
must be sensitive to the staggered field. This means that 
it can only contain bulk solutions which are substantial 
superpositions of k and fc + tt, i.e. b^/aj: ^ 1. From the 
expression, Eq. for 6g/ag we find = ^{a/b — b/a)m 

which implies that k is constrained to the volume |e^| ^ 
m. 




FIG. 2: Spectral weight of various charge potentials, (a) is a 
stripe along the y direction with m = 2,t' = ~.l,t" — .1, (b) 
is a diagonal stripe along x + y direction with m = 0.5, t' = 
— .l,t" = and (c) is a size 10 x 10 bubble with m = l,t' = 
-.l,t" = 0. 

In Fig. 121 is shown the full spectral weight of the three 
types of potentials displayed in Fig. ^ confirming our 
analytic result. The broader distribution in k-kspace for 
larger m is consistent with a shorter localization length 
in real-space. Note also that for the diagonal stripe. Fig. 
13), the spectral weight is concentrated to those segments 
of the diamond which have momenta orthogonal to the 
stripe direction, a feature which will be important when 
analyzing diagonal stripes in the next section. 



III. DIAGONAL STRIPES 

One of the most interesting recent stripe related find- 
ings is the diagonal stripes seen by quasi-elastic neu- 
tron scattering in very underdoped insulating phase of 
LSCO.^^ As discussed in Sec. [Jthere is a crossover region 
where both diagonal and bond-aligned stripes appear to 
coexist and which coincides with the rapid drop of the 
superconducting transition temperatures with decreasing 
dopingfi^ pointing to a strong connection between stripes 
and superconductivity. On the other hand, we also noted 
that other properties, such as the nodal Fermi velocity 
and the hole mobility does not show any dramatic change 
across the bond-aligned to diagonal transition. 



What we can contribute to this discussion within our 
model is a comparison between the expected distribu- 
tion of spectral weight between bond-aligned and diago- 
nal stripes. It was already shown in Ref. 11 an dT? that a 
disordered bond-aligned stripe array can well reproduce 
the qualitative features of the near optimally doped sam- 
ples. It is natural to do a similar analysis for the lightly 
doped samples a; = 3 — 5% studied by ARPES in Ref. 
lid Experimentally it is found that the spectral weight 
in the anti-nodal, (7r,0), regions of the BZ which is most 
prominent for the under and optimally doped samples is 
gapped away from the Fermi energy, and instead the low- 
energy spectral weight consists of disjoint arcs of "Fermi 
surface" centered near the nodal, (7r/2,7r/2), regions. At 
a first glance it is very tempting to identify these features 
with that shown in Fig. for a diagonal stripe. How- 
ever, although it is probably not completely ruled out 
that this simple picture is correct a closer study of the 
diagonal stripe reveals a serious problem. The problem 
is that for t' < 0, a. less then half-filled, i.e. hole doped, 
diagonal stripe will have very little spectral weight in the 
nodal region, but only near the anti-nodal regions. 

In order to understand this statement we can consider 
the dispersion of the diagonal and bond-aligned stripe as 
a function of the conserved momentum along the stripe. 
We will restrict ourselves to the limit to — > cxd and nu- 
merically confirm the qualitative correctness for to = 1. 
In this limit it is trivial to solve for the spectrum on the 
stripe because the problem reduces to a one dimensional 
tight-binding chain. For the bond-aligned stripe there is 
the nearest neighbor hopping t while t" acts as a next- 
nearest neighbor hopping, resulting in a dispersion 

£^„i(A:||) '"=°° -2icos(fc||) - 2t"cos(2A:||) . (8) 

For the diagonal stripe the next-nearest neighbor hopping 
on the 2D lattice t' acts as a nearest neighbor hopping 
on the chain, giving a dispersion 

ediag(A:||) "=°° -2t'cos(A:||) . (9) 

In Fig. 131 are shown numerically calculated dispersions for 
TO = 1 which agrees qualitatively with the large m limit. 
Note that for the bond-aligned stripe there is a folding of 
the BZ due to the antiferromagnetic scattering along the 
stripe direction, which is absent for the diagonal stripe. 
The conclusion we want to draw from these dispersions 
is that for the diagonal stripe at any finite hole doping 
the momenta of filled states which may contribute to the 
low energy spectral weight are confined to > 7r/2. 
From the distribution of spectral weight for a diagonal 
stripe. Fig. |3), we find that this implies that there is 
very little spectral weight in the nodal regions. We can 
easily convince ourselves that this distribution of spectral 
weight of a diagonal stripe is a general consequence of the 
fact that a state with some momentum fcy will, as shown 
in Sec III Al have its spectral weight concentrated to the 
intersection of the line k^ + ky = fc|| with the cos(fc2,) -I- 
cos(fcj,) = diamond. This implies that for |fc||| = tt the 
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weight will spread out over the whole intersection which 
is a line, while for |fc|| | < tt the weight will be concentrated 
to the two points of intersection. 
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FIG. 3: Dispersions of in-gap states for m — t — 1, t' = —.1, 
and t" = for bond-aligned stripe (a) along y with fcy — ky 
and diagonal stripe (b) along x + y with fcy — kx + ky. Line 
marked "0%" indicates no doped holes (half- filling). 

These results for the distribution of spectral weight of 
a hole doped diagonal stripe may be contrasted with that 
of a bond-aligned stripe. Here we find fromlSJ^ comple- 
mented with Fig. that for finite hole-doping the spec- 
tral weight may be spread out over the whole diamond. 
The details, of course, depending on the doping, the pa- 
rameters used and on the density of stripes. The latter 
being particularly important in that for a short inter- 
stripe distance the stripe states will overlap and form 
bands of momenta transverse to stripes (see Fig. I14|) . 

One might, in an effective model such as this, attempt 
to fit the experimental data by taking t' > which would 
allow for spectral weight concentrated in the nodal re- 
gion. However the ARPES data for the lightly doped 
samples indicate that most of the additional weight in- 
troduced with doping is in fact in the anti-nodal region 
only that it is gapped away from the Fermi surface. If we 
take i' > we would vacate the anti-nodal states and we 
would not be able to reproduce this qualitative feature. 
Related to this there is a more quantitative problem for a 
diagonal stripe as contrasted with the ARPES data which 
should be more general than our model, namely that the 
band width of a purely diagonal stripe is expected to be 
proportional to t'. Values of t' in the literature is less 
than 0.1 eV, implying a band width Wdiagonai < -2 eV, 
whereas the band width from the in-gap states seen in 
ARPES can be estimated at Win-gap ^ 1 eV which looks 
more consistent with the band width Wcoi 2f of bond- 
aligned stripes. To summarize, we find that pure diagonal 



stripes are not consistent with the ARPES data of lightly 
doped LSCO. 

A. Staircase Stripes 

Given the difficulties with matching the model using a 
diagonal stripe configuration to the ARPES data it is in- 
stead tempting to look at bond- aligned stripes. Now, we 
know from neutron scattering that bond-aligned stripes 
are not seen in these very lightly doped materials, but 
only diagonal stripes. This led us to investigating the 
properties of stripes which are locally bond-aligned but 
globally diagonal. A natural and most simple candidate 
for such a construction is a "staircase" stripe. We can 
define a staircase stripe along the x + y direction but let- 
ting it run alternately along the x and y directions with 
some step length I. For an ordered array of such stair- 
case stripes we also introduce a stripe distance d defined 
according to Fig. 0| In the case of anti-phase stripes the 
potential has the symmetries 

V{x + l,y + l) = V{x,y), 
V{x + d,y-d) - -V{x,y), (10) 

which also gives the primitive cell as indicated in the fig- 
ure. The same symmetries hold true for the full potential 
of Eq. n which is simply multiplied by a factor (— 1)^+^ 
to account for the staggered field. 
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FIG. 4: Graphical representation of a staircase stripe defined 
by the step length I and the stripe distance d. The rectangle 
indicates a primitive cell. The potential V(x,y) is given by 
V=0 for darker gray, V=l for light gray, and V=-l for white. 



We will return to magnetic structure factors of such 
staircase stripes below, but it is easy to see that as long 
as I < d the main magnetic diffraction peaks of such 
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a staircase stripe are equivalent to an array of purely 
diagonal stripes with the inter stripe distance 2d along 
the X and y direction. 

We turn now to the distribution of spectral weight 
of staircase stripes. As an example we look at at a 
system with / — 8 and d = 8, were d is chosen such 
that the magnetic structure factor has main peaks at 
(tt ± (5/V2,7r =F 5/^/2) with 5/^/2 = 1/32 which corre- 
sponds roughly to the 5 w 1/25 seen in Neutron scat- 
tering at 4% doping. We diagonalize this system numer- 
ically to find the single particle eigenstates ipaik) with 
energies Ea in terms of which we calculate the single 
particle spectral function 



■UJ 



(11) 



where 6{Ea — uj) is the Kronecker delta, and the local 
density of states 



(12) 



with ipair) = 1/LxLy J^j^ e ipa{k) being the eigenfunc- 
tions in real space. In Fig. |5l is shown spectral weight 
distribution in fc-space and real-space when integrated 
over an energy window Aw — .2 around the Fermi en- 
ergy at 4% doping, i.e. calculating 



and 



pEf+Auj/2 

I{k) = / A{k,oj)duj, 

J_Ejr-Aw/2 



R{r) = I R{r, Lo)dLo . 



(13) 



(14) 



We find that the low-energy spectral weight is concen- 
trated near the nodal region. It is highly anisotropic with 
most of the spectral weight parallel to the overall stripe 
direction along x + y in contrast to the pure diagonal 
stripe shown in Fig. ^p. We have tried making the step 
length / shorter, this gives results closer to the pure diag- 
onal stripe with most of the spectral weight in the anti- 
nodal region. In Fig. Elthe results are symmetrized with 
respect to the stripe direction, so that (tt, tt) = (tt, — tt) 
etc. Here the disjoint features merge into a single piece of 
"Fermi arc" in each quadrant of the BZ, similar to what 
IS seen m ARPES. Note that in Fig. ^ it is not the stripe 
potential V{x,y) which is plotted, but the amplitude in 
real space of the low energy states. Not surprisingly these 
follow the potential quite closely, but there is some leak- 
age of spectral weight into the antiferromagnetic regions 
which appears to smooth the kinks and make the stripes 
more diagonal. An indication that the staircase stripes 
are just a caricature with the real stripes probably being 
smoother, but nevertheless locally closer to bond-aligned 
than to diagonal. 

We can gain a better understanding of these results for 
the low-energy spectral weight by studying the spectral 





(a) 




r 








* k. 



FIG. 5: Ordered array of staircase stripes with / — 8 and 
d — 8, using t — m — 1, t' = —.1 and t" — 0. Momentum 
space (a) and real space (b) spectral weight integrated over 
an energy window of 0.2 around the Fermi energy at 4% hole 
doping. The full system size is 256 x 256 while in (b) is shown 
a 100 X 100 section. 
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FIG. 6: Same as Fig. I^but symmetrized with respect to the 
stripe orientation and in the first quadrant of the BZ. Shown 
is also the intensity map which has a linear scale. 



hmction, A{k,Lu), over a broader energy window along 
the high- symmetry directions, as shown in Fig. [3 The 
qualitative features of the spectral weight deriving from 
the staircase stripes can be directly linked to the prop- 
erties of ordinary bond-aligned stripes. The high energy 
(away from Ep) spectral weight is concentrated around 
the anti-nodal, (tt, 0), region whereas the low energy spec- 
tral weight close to the Fermi energy is focused to the 
nodal region around (7r/2,7r/2). This is what we would 
find from a lightly hole doped bond-aligned stripe with 
the dispersion in Fig. Oi, given the fact that the spectral 
weight is concentrated to the intersection of /c[| = or 
ky with the BZ diamond. 



1. Magnetic Structure factors of Staircase Stripes 

In the previous section we found a qualitative agree- 
ment of the spectral distribution of staircase stripes with 
that seen experimentally in lightly doped LSCO. Here 
we had to restrict ourselves to staircase stripes with step 
length / less than the stripe distance d as defined in Fig. 
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FIG. 7: Band structure of system with t — m = 1, t' — — .1 
and t" = for (a) no stripes and (b) with staircase stripes as 
in Fig. 13 UHB and LHB indicate upper and lower Hubbard 
bands respectively and Ef is the Fermi energy at 4% hole 
doping. The spectral weight is indicated by the intensity, but 
on a non-linear scale which exaggerates low-intensity features, 
(b) is symmetrized with respect to the stripe direction. 

01 in order to have a magnetic structure factor which cor- 
responds to diagonal stripes. A natural extension is to 
study also stairs with I > d. Certainly, in limit I 3> d 
we expect very little influence from the kinks and the 
system will become equivalent to that of bond- aligned 
stripes, both for the spectral distribution and structure 
factors. We found that the distribution of spectral weight 
could be rationalized in terms of bond-aligned stripes al- 
ready for the system with I = dso that for ^ > d we would 
not expect any qualitative difference between the spectral 
distribution of staircase stripes with that of bond-aligned 
stripes studied in the earlier workiii*i^ What is more in- 
teresting is to study the magnetic structure factor, which 
is sensitive to the global properties of the system. 

Physically, the relevant entity is the magnetic structure 
factor which is the amplitude |S'^((7)p of of the Fourier 
transformed spin density which at zero temperature reads 

S'iq) = J2 e''-V|Va.(r)|'e(i?j. - E^) . (15) 

r,a,cT 

However, for a large system, such as for disordered 
stripes, this is difficult to calculate because of the need 
to diagonalize the system. Much simpler to find is the 
amplitude squared of the Fourier transform of the stripy 
potential {—l)^~^yV{x,y). We have checked for several 
ordered stripe arrays that close to half-filling the two 
calculations give very similar results. This is expected 



because the difference is roughly the amplitude on the 
stripes which is zero for the potential but slightly differ- 
ent from zero for the actual spin density because of the 
occupied stripe states. 

Using the symmetry properties of a staircase stripe, 
Eq. 1101 wc find that Vikx^ky) can only have non-zero 
components for 

k,j,^ky = —N, N e Integer (16) 

k^-ky = — iV', TV' e odd Integer . (17) 

In Fig. |S1 is shown the structure factors of three dif- 
ferent stripe realizations, I = d, I = 2d and I = 5d. We 
find three very distinct diffraction patterns, were the first 
corresponds to diagonal stripes, the second looks like di- 
agonal stripes coexisting with both vertical and horizon- 
tal bond- aligned stripes, and the last looks like vertical 
and horizontal stripes but with an orientation which de- 
viates slightly from bond-aligned. We also show diffrac- 
tion patterns of samplings over disordered stripe config- 
urations were we find that in general only the primary 
peaks survive with secondary, lower intensity peaks, get- 
ting washed out, although this does not happen for the 
second configuration with coexisting bond-aligned and 
diagonal peaks. We have also studied stripe arrays with 
larger unit cells where I is not an integer factor of d. For 
disordered realizations it appears that these roughly fall 
into one of the three characteristic regimes, with only a 
quite narrow window of Z « 2d showing both diagonal 
and bond-aligned peaks of appreciable amplitude. 

These results for the diffraction pattern of staircase 
stripes appear qualitatively very similar to what is seen 
in neutron diffraction experiments on LSCO in the low- 
temperature orthorhombic (LTO) phase.— In the non- 
superconducting phase at very low doping {x < 6%) is 
found peaks consistent with diagonal stripes, whereas at 
higher doping (a; « 10 — 13%) are found a pattern con- 
sistent with stripes that are close to bond-aligned but 
shifted by a few degrees from the tetragonal axes. More 
recently was found that very close to the insulator to su- 
perconducting transition (x ~ 6%) both patterns coex- 
ist. This may suggest that the stripes in the LTO phase 
of LSCO are always of the staircase type and that the 
transition between bond- aligned and diagonal stripes is a 
crossover from a regime where I > 2d to one where I < 2d. 
Particularly interesting would be to have more detailed 
experimental data from neutron scattering on the un- 
der doped superconducting regime x ~ 10%. From the 
staircase stripe scenario we would expect a larger angle 
of deviation from the tetragonal axes than for the sam- 
ples with higher doping, possibly together with secondary 
peaks in the diagonal direction or on the tetragonal axes 
as shown in Fig. |21 

We should note that the correspondence between ex- 
periment and our results for staircase stripes is not per- 
fect. In particular there appears to be a discrepancy 
with the relation between the "incommensurabilty" S, 
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FIG. 8: Magnetic structure factors of the three different 
regimes of staircase stripes discussed in text. The left column 
shows the unit cells which in (a) is 80 x 80 with step length 
I = 20 and stripe distance d = 20, in (d) is 40 x 40 with 
/ = 20, d = 10 and in (g) is 80 x 80 with / = 20, d = 4 . Black 
are the stripes where the potential is zero, gray and white 
are the two vr-shifted domains of the antiferromagnetic order, 
and the rectangles are primitive cells. The middle column, 
(b), (e) and (h), shows the corresponding structure factors in 
reciprocal lattice units centered around (1/2,1/2) = (7r,7r). 
The right column, (c), (f) and (i), give the structure factors 
of a sum over ten size 80 x 80 disordered configurations of the 
corresponding ordered states to the left, where the length of 
the legs and the distance between stripes is allowed to vary 
with a flat random distribution without allowing stripes to 
come closer than one site separation. The width of the peaks 
of the right column is simply related to the fraction of the 
BZ viewed, given that in all calculations the same correlation 
length is used. 

defined as the shift of a peak from (vr, tt), of the diagonal 
and bond-aligned components. For the staircase stripes 
we find (Jdiag = ^coi/ V^, whereas experimentally is seen 
Sdiag ~ i^coi- The former relation follows directly from 
the symmetry of the stripy potential, Eq. 1171 and conse- 
quently is not sensitive to small changes of the potential 
such as for instance smoothing of the kinks. 

IV. IN-PHASE STRIPES 

We now turn to a study of the spectral weight distribu- 
tion of in-phase stripes. As discussed in Sec. Qlthere is in- 
direct evidence for an inhomogeneous charge distribution 
in the electron doped PLCCOii^. In fact, the possibil- 
ity of in-phase stripes had already been suggested from 
theory for hole doped stripes at lower hole densitieSi.2i 
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FIG. 9; Magnetic structure factors of disordered realizations 
of staircase stripes with I = 20 and (a) d = 5 and (b) d = 6 
which could roughly correspond to 10% and 8% doping re- 
spectively. 



The experimental hallmark of anti-phase stripes is the 
incommensurate magnetism detected by neutron diffrac- 
tion, where the weight is shifted away from the antiferro- 
magnetic ordering vector Q = (tt, tt). A system with or- 
dered in-phase stripes on the other hand will have a main 
diffraction peak at the antiferromagnetic ordering vector 
with satellites at positions shifted by integer multiples of 
27r/(i, where d is the distance between stripes. However, 
if the stripes are not static and ordered, but fluctuating 
or disordered, the satellite peaks will easily be washed out 
and the weight absorbed into the AF peak. In Fig. 1101 is 
shown the amplitude squared of the Fourier transform of 
the stripe potential {—l)^'^^V{x,y) along the transverse 
direction for ordered and disordered anti-phase and in- 
phase stripes and the for the charge density represented 
by For the charge order there is no difference 

between anti-phase and in-phase stripes. Nevertheless, 
charge order is more difficult to detect even for systems 
with static order. The reason for this is that neutrons do 
not couple directly to the charge order and that the su- 
perlattice peaks related to charge stripe order arise from 
modulations of the uniform charge density, Q = (0,0), 
which dominates the structure factor. As shown in Fig. 
[m these peaks are also readily destroyed by disorder. 
(In fact, the charge structure factor as estimated here is 
equivalent to the spin structure factor of in-phase stripes 
but shifted by tt.) Direct signatures of in-phase stripes 
from neutron or x-ray scattering will thus be much more 
difflcult to find, most likely requiring a stripe ordered 
material. 

Nevertheless, given the indirect indications for the ex- 
istence of in-phase stripes in electron doped cuprates 
it may be interesting to study the implications of such 
structures to the distribution of spectral weight. In Fig. 
[TTI is shown the integrated spectral weight close to the 
Fermi energy for disordered in-phase stripes with inter- 
stripe distances 9-15 with mean 12 and 1-7 with mean 4 
for a system with t = m — 1, t' — —.2 and t" ^ .1. The 
parameters t' and t" are chosen such as to roughly repro- 
duce the band structure of the undoped system seen in 
ARPES on the electron doped NCCO^ with the small- 
est spectral gap at (7r/2,7r/2) as shown in Fig. I12h ,. The 
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FIG. 10: Structure factors (in arbitrary units) of spin and 
charge density for anti-phase and in-phase stripes as a func- 
tion of momentum transverse to the stripes in r.l.u.. The left 
column shows ordered arrays with stripe distance d = 4 and 
the right disordered with a flat distribution d — 2 — 6. (a) 
and (b) shows the spin order for anti-phase stripes, (c) and 
(d) the spin order for in-phase stripes, and (e) and (f) the 
charge order which is independent of the type of stripe. In all 
figures is used a finite correlation length of 40 sites and the 
disordered systems are averaged over IQQQ samples. 
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FIG. 11: Spectral weight close to Ef (window .2) of dis- 
ordered in-phase stripes with t = m — 1, t' — —0.2, and 
t" — 0.1. Mean stripe distance in (a) is 12 with Ep at 4% 
doping and in (b) mean distance is 4 with Ef at 10% dop- 
ing. The results are symmetrized with respect to the stripe 
direction. The system size is 320 x 320 . 



density of stripes is assumed to be roughly as in the hole 
doped materials with stripe spacing d given by the dop- 
ing n according to n « l/2oJ, which for electron doping 
corresponds to 3/4 filled stripes. 

Clearly these results reproduce the experimental Fermi 
surface quite well, at light doping there are patches of 
spectral weight around (tt, 0) while at higher doping the 
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FIG. 12: Band structure of system with t — m — 1, t' = —0.2, 
and t" = 0.1 and in-phase stripes as in Fig. 1121 (a) is no 
stripes, (b) mean period 12 stripes and (c) mean period 4 
stripes. 



weight starts to look more like a full Fermi surface closed 
around (tt, tt) but with weight missing at "hot-spots" 
where the putative Fermi surface cuts the BZ diago- 
nal (diamond). Similar results have been reproduced by 
mean-field theory of the t — t' — t" — U Hubbard model 
by allowing for a doping dependent U .^^ This is in fact 
the model we consider but without the stripes and with a 
magnitude of the staggered field m which depends on U. 
Here, we keep the parameters fixed but vary the stripe 
spacing as a function of doping. There are however very 
distinct differences between the implications of the two 
scenarios. In the stripe model there are mid gap states, 
not present for the model with a uniform staggered field. 
This is particularly visible for the lightly doped system, 
Fig. 112b . where along the F = (0,0) to (tTjTt) direc- 
tion there is spectral weight between Ep and the lower 
Hubbard band. Precisely such a feature is seen in the 
ARPES data for the 4% doped sample. Secondly, the 
low-energy states of the stripe model are dynamically 
one-dimensional, i.e. they are localized transverse to 
a stripe but have a well defined momentum along the 
stripe, as shown in Fig. 1131 Thus, one of the main 
conclusions is that the distribution of low-energy spec- 
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tral weight, the "Fermi surface", may be indistinguish- 
able from a homogeneous two-dimensional system even 
though the states are dynamically one- dimensional. 



.2 




(a) 


. 1 






k 


... ...xk 





. 02 




100 200 300 



7T/2 



37r/2 



FIG. 13: Example of a single state of system in Fig. I12t . 
with energy E ~ and momentum along stripe direction ky — 
it/2, showing localization in the transverse stripe direction. In 
(a) is shown the real space amplitude \ip{x,y)\'^ as a function 
of transverse direction x for arbitrary y and in (b) the k-space 
amplitude We can estimate a localization length 

^ « 2 from the real space peak amplitude. 



Let us now look at how we can understand this evolu- 
tion of spectral weight with doping and corresponding in- 
creasing stripe density from the properties of a single in- 
phase stripe. Fig. 1141 shows the spectra of systems with 
sparse, mean distance 12, and dense, mean distance 4, 
stripes. For the sparse stripes, Fig. [Hk . the stripe states 
are clearly visible as the isolated in-gap bands. Because 
there is no broadening of the stripe states the stripes are 
clearly sufficiently far apart as to be effectively indepen- 
dent. We find that in-phase stripes at moderate m have 
a band structure which is distinctly different from that of 
anti-phase stripes as shown in Fig. with the former 
staying close to the upper and lower Hubbard bands. (We 
have checked that in the limit m — s- oo, where in-phase 
and anti-phase stripes are equivalent, the dispersion Eq. 
ISjis correctly reproduced also for in-phase stripes.) As 
the stripe density is increased the stripe states will over- 
lap and form bands transverse to the stripes. Fig. 114b . 
with states from the lower branch crossing the Fermi en- 
ergy. Because of disorder, these states may nevertheless 
be strongly localized as shown in Fig. 1131 

An important point about in-phase stripes which is 
demonstrated here is that the chemical potential may 
move as a function of doping and is not necessarily 
pinned within the gap. This is particularly evident in 
comparing the 0% and 4% samples. Fig. [T^ and b, 
where in the former the chemical potential is in the gap 
(—0.6 < /i < —0.2) whereas in the latter it has moved 
up to cut the upper Hubbard band {fi ~ 0). The motion 
of the chemical potential is a necessary consequence of 
the dispersion of in-phase stripes as shown in Fig. I14h ,. 
with the stripe states "hugging" the upper and lower 
Hubbard bands. This is in sharp contrast to anti-phase 
stripes (Fig. [7|and|3J), where the "stripy" spectral weight 
is "mid-gap" implying that the chemical potential may 
stay fixed with doping. Motion of the chemical potential 
with doping is thus not necessarily an indication of the 
absence of stripes P'^'i^^ 
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FIG. 14; Spectrum for \E\ < 2 of in-phase stripe systems as 
m Fig. iniwith mean stripe period 12 in (a) and 4 in (b) as a 
function of momentum k\\ — ky parallel to the stripe direction 
y. The black curves indicate the bounds of the upper and 
lower Hubbard bands for the corresponding system without 
stripes. 



A. Bubbles 

Because of the lack of direct diffraction evidence for 
stripes one may be free to speculate on other forms of 
charge order in these materials. The microscopic mo- 
tivation for stripe formation is the tendency of the an- 
tiferromagnet to expel extra charge which will disrupt 
the local antiferromagnetism. The formation of stripes is 
then a compromise between the minimizing the magnetic 
energy by concentrating the holes, kinetic energy by al- 
lowing holes to delocahze along the stripes, and possibly 
the charging energy due to long range Coulomb repulsion 
by not allowing a macroscopic charge inhomogeneity.^^ 

However, if the magnetism is relatively stronger than 
the kinetic energy contribution the system may prefer to 
keep the extra charge in zero-dimensional "bubbles" (see 
Fig. ^) in order to minimize the disruption of the lo- 
cal magnetic order i2i Roughly speaking, for a completely 
filled stripe with two electrons per site on a one-site wide 
stripe there are three bad bonds where the spin exchange 
is destroyed per stripe site, while for a bubble there are 
only two bad bonds per site. (The sites on the perime- 
ter of the bubble have three bad bonds). In addition, 
because of charging energy, the bubbles would have to 
be limited to a microscopic size. Clearly, the putative 
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stripes or bubbles are not completely filled in the elec- 
tron doped materials, due presumably to kinetic energy 
considerations, as this would imply an insulating system. 
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FIG. 15: Spectral weight within energy window .2 of Ef of 
ordered array of bubbles with t — m = 1, t' = —0.2, and 
t" = 0.1. In (a) the unit cell is 12 x 12 with the bubble size 
4x4 and Ef is at 4% electron doping. In (b) the unit cell is 
20 X 20 with the bubble size 10 x 10 and Ef at 10% doping. 



For simplicity we look at square bubbles and because 
this problem is fully two-dimensional with a large unit 
cell we only consider ordered systems. For the previ- 
ous study of in-phase stripes we chose stripe spacings 
and fillings which were based on the corresponding val- 
ues for anti-phase stripes in the hole doped systems. For 
the bubble phase we have even less to guide us on how 
to choose the size and period of the bubbles. However, 
assuming that the putative bubble formation is due to 
competition between the magnetic energy and the charg- 
ing energy we can get a simple estimate of the variation 
of bubble size with doping, ignoring in this completely 
the kinetic energy. It is easy to see that for a classical 
antiferromagnet we decrease the number of bad bonds 
which are not connecting antiferromagnetically aligned 
spins by two per doped electron by forming bubbles in- 
stead of a homogeneous (Wigner crystal) distribution of 
doped charge. This gives an energy gain Eaf ^ —JL'^, 
where J is the AF exchange which is typically of the 
order of 0.1 eV and L is the linear dimension of the bub- 
ble. The charging energy due to moving electrons 
together from a closest distance 1/n to a closest distance 
1 is given by Eq ~ QL^{l—n) assuming a Coulomb inter- 
action V = Q/r with r in units of the lattice constants. 

2 

Here Q = -r-^ — « .3 eV assuming a lattice constant of 

a ~ sAand a dielectric constant e = 10.^^- Minimizing 
the total energy with respect to the bubble size L gives 
so that for small doping, n, we find the bub- 



(l-n) 



ble size proportional to the doping. The prcfactor J/Q 
seems somewhat on the small side but does not rule out 
bubbles as a viable scenario. Clearly large J, or "strong" 
antiferromagnetism, will favor the bubble scenario. This 
being the reason why we suggest it for the electron doped 
materials which in general show antiferromagnetism over 
a much wider doping range than do hole doped materials. 

For the calculations presented here we choose a bubble 
size L = lOOn and let the distance between bubbles vary 
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FIG. 16: Band structure of system in Fig. 1151 The band 
structure without bubbles is the same as in Fig. 112b . 



with doping such as to correspond to roughly 3/4 filled 
bubbles. In Fig. El is shown the the integrated spec- 
tral weight near the Fermi energy for two realizations, 
at 4% and 10% doping, and in Fig ^1 the corresponding 
band structure. The results are qualitatively similar to 
what was found for in-phase stripes. The evolution of the 
weight near (tt, 0) does not depend crucially on the bub- 
bles but can be understood from increasing filling of the 
upper Hubbard band. The nodal weight near (7r/2, tt/2) 
on the other is clearly an effect of the additional states 
deriving from the bubbles. 



V. CONCLUSIONS 

We have investigated the distribution of electronic 
spectral weight in various charge ordered antiferromag- 
net s. 

We find that the spectral weight of states localized 
on a stripe or other charge structure is centered on and 
spread out over the Brilloinn zone diagonals {cos{kx) + 
cos(fcy) = 0). The distribution of spectral weight close to 
the Fermi energy may thus look fully two dimensional and 
practically indistinguishable from a homogeneous system 
even though the low-energy states are dynamically one- 
or even zero-dimensional, localized on stripes or bubbles. 
On the other hand the stripe states will in general lie 
within the energy gap of the undoped system making the 
appearance of "in-gap" states an expected consequence 
of an inhomogeneous charge distribution. 

We find that pure diagonal stripes cannot reproduce 
the distribution of spectral weight found in ARPES on 
the very underdoped LSCO in the "diagonal stripe" 
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phase. Instead we introduce "staircase" stripes which 
are locahy vertical or horizontal but globally diagonal, 
in terms of which the qualitative features of the ARPES 
data is readily reproduced. Calculating the structure fac- 
tors of such staircase stripes we find that these evolve 
with doping and corresponding stripe density in a way 
which is very similar to the neutron scattering data in 
the LTO phase of LSCO over the whole doping range 
from very light to optimal. The results suggest that 
that the horizontal to diagonal stripe transition may be 
a crossover between a regime where the typical length of 
straight (horizontal or vertical) stripe segments is longer 
than the inter-stripe spacing to one where it is shorter 
and that the stripes are always locally bond-aligned in 
LSCO. 

We find that in-phase stripes can qualitatively repro- 
duce the ARPES data end evolution with doping in the 
electron-doped cuprate NCCO. Particularly revealing is 



the spectral weight in the nodal region, near (7r/2,7r/2), 
where at low doping there is "stripy" spectral weight in 
the gap which at higher doping broadens as a result of 
increased stripe density and crosses the Fermi energy. 
For in-phase stripes, in contrast to anti-phase stripes, we 
find that the in-gap states lie close to the upper and lower 
Hubbard bands implying that the chemical potential is 
likely to move with doping. We also consider bubble 
structures and find that these produce similar results for 
the distribution of spectral weight to that of in-phase 
stripes by allowing the bubble size to grow with doping. 
We argue that these may be an alternative to stripes in 
the electron doped materials where there is a broad dop- 
ing range with antiferromagnetic order. 

We want to thank S.A. Kivelson for valuable discus- 
sions. This work was supported by the Swedish Research 
Council. 
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